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I'lntroduction 

An  impottant  enoneoing  proUem  is  the  detenxdnatkm  the  electromagnetic  fields 
in  miaowave  systems,  m  example  &q)ered  waveguides,  boms,  scattercrs,  clo^  cavities, 
and  (^)ai  resonators.  We  heg^consider  the  case  oi  axisymetric  transvase  dectric  modes. 
Such  problems  for  monochrooiuic  radiation  can  be  reduced  to  considentdon  of  an  ellipdc 
partial  ifiEFerential  equation  amilar  to  die  Hdmhdtz  equatioiPTMedKxls  have  been 
deyelqped  for  the  direct  numerical  soludon  of  the  partial  differential  equatjon^TVariatirMal 
princyles  have  been  taed  to  optinMdlydeterajineiCTBMdinate  values  Adjects  of  interest 
like  reflection  and  transmission  coefi^entsT^  alternative  approadi  is  die  reduction  of 
the  problem  ^consideration  of  an  int^ral  equation  defined  on  die  metallic  walls  defining 
the  object boundary  integral  methodVThjBse  have  been  sdved  for  the  case  of  scalar 
fidds  described  by  the  Helmholtz  equatio^vlbe  boundary  integral  equation  method  is 
feasible  when  die  Greens  function  is  known  in  a  cornputationaly  convenient  form,  and  is 
very  ofien  much  more  cooqiotatiooaly  effident  than  its  oon^titors,  p^cularly  when  the 
geometry  is  conqilex.  Ws  here  dcrwiBCthe  dieory  and  effective  numerical  implementatirxi 
of  such  a  boundary  integral  equationi^roach  for  die  case  of  an  axisymetric  transverse 
dectric  electromagnetic  fidd.  The  tech^ue  is  readily  generalizable  to  arbitrary 
axisymmetric  fields.  A  »  ;:>  - 

y 

2  -  Mathematical  and  Phydcal  Prdiminarics 

Considera  transverse  electric  field  which  in  cylindrical  coordinates  p,^,z  can  be 
expressed  as 


r  ‘ 


E»(V<^)Re[e-*®t‘P(p,z)]  (1) 

Notethat  V^se^p.  The  Maxwell  equaticHis  read  for  this  case 

VxE  « icopH  (2) 

VxH^-ioeE  (3) 


where  the  didectric  constant  e  and  pennedaJity  p  will  be  taken  to  have  dieir  vacuum 
values.  The  situation  where  they  ate  piecewise  constant,  appropriate  for  example  for  a 
waveguide  with  didectric  Uy^  can  readily  be  treated  ^  a  generalization  of  die  methods 
to  be  presented.  It  fdlows  dire^y  firona  (1)  and  (2)  that 

n  j  i 

eop  p«op 

If  one  forms  die  scalar  product  of  (3)  with  ,  on  using  (4)  to  eliminate  H  and  (2)  to 
express  E,  there  results  after  some  algebra 
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or  on  writing  dungs  out  in  cylindrical  cooRUnates 


The  tinoB  avenged  Poynting  vector 


r«i 


4agip2 


and  it  is  readily  established  fiom  (5)  that 


(5) 


(6) 


(7) 


vr  =  0  (8) 

Eq.  (7)  can  ctuiveniendy  be  enq>loyed  in  testing  numerical  results. 

In  vacuum  xl/c^,  vritere  c  is  the  speed  of  light  Define  ksoVc.  We  shall  require  for 

what  follows  the  ou^oing  Greens  function  anociated  with  (S).  It  is,  as  shown  in 
Appendix  A, 


K(p,p',z-z')  =  pp’ 


e -»(♦-♦’) 


where,  since  r  «  pep-t-zcz  and  r'=  p'CpM-z'ez » the  separation 


(9) 


R  =  Ir - r'l  =  [p2  +  p'2  -  2pp’cos(4>  -  ♦’)  +  (z  -  z’)^]*^ 
.  As  is  also  shown  in  qtpendix  A 


L{  K) 


8(p-p*)S(z~z’) 

(pp*)i/2 


(10) 


(11) 


Note  that  K  is  invariant  under  interchange  of  primed  and  unpiimed  arguments.  Thus 


1  a  /I  8Kn  d  /  I  aKx  kZK  S(p  -  p’)S(z  -  zQ 

p’dp'Vp' 3p'J  p*2  (pp’)i^ 


L’{K} 


(12) 
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The  Greens  function  can  be  used  to  reduce  the  problem  of  determining 'F  to  the 
stdution  ol  an  integral  equation  on  die  perfecdy  conducting  metallic  surfaces  assumed  to 
bound  the  system.  This  serves  to  replan  die  partial  differential  equation  (5),  definedcma 
two  dinaentiooal  domain,  by  a  proUem  defined  on  a  one  diinensioDal  domain.  There  are, 
indeed,  a  number  of  such  intq^  equations.  We  shall  deal  with  one  such  which  is  the 
second  kind  and  provides  a  formulation  which  is  nmnerically  stable. 

Suppose  that  die  domain  interest  is  a  volume  V  bounded  by  a  surfece  S  obtained  by 

rotating  around  die  z  axis  a  curve  r  defined  in  a  plane  ^iaconst  Let  sdenote  arc  length 
along  r ,  and  suppose  diat  F  is  defined  by  the  parametric  equatims  psp(s)  and  zsz(s) . 
Then  dss[(tip)^tiz)^^.  The  volume  element  is  (Prspdpd^idz.  The  vector  donent  ctf 
surface  is  d^rsnrdsdi^ ,  where  n>mp(p,z)Cp-Hiz(p,z)ez  is  die  unit  normal  pointing  out  of  V. 
We  seek  an  axisymmetric  function  a(p,z)  defined  on  A  such  that 

(13) 

or  equivalently 

'P(p*,z')  =  Jds  (14) 

r 

PCX'  any  point  p',z'  interitx  to  V  and  not  on  Fit  follows  from  (1 1)  and  (13)  that 
L'('F(p',z’)}sO.  But  on  F,  as  shown  in  appendix  A,  because  of  the  singular  nature  of  K 
as  r-^r*,  eq.  (13)  reduces  to 

♦(s')  =  o(s')  +  jF(s',s)o(s)ds 
F 

where  the  inhomogeneous  term 


is  presumed  to  be  given.  The  kemd 


F(s' s)  =  ii(p  z)-V  ^ 

(S.S)  (P.  )  1^  Jp«p(s)  z«z(s)  p'»p(s')  z’»z(s') 


(17) 


deveh^  a  ample  pole  as  s-m',  whence  the  integral  noust  be  interpreted  as  a  Cauchy 
{mncipal  value.  The  factor 


•srr 


(18) 


Eq.  (IS)  is  the  desired  integral  equation  of  die  second  kind  for  o. 

For  the  case  of  horns  and  scatterers  one  iscoocented  with  die  field  far  fiom  the  object 
This  is  readily  determined  fiom  the  integral  lepresentadmi  (13).  To  this  end  note  that  as 

r’=<p’2+z'2)l/2  — >  oo,  one  has  Ir-r*!  ~  f  -  r*  (r’A')  +•••,  ^ence  to  lowest  significant 
(miff  in  T/f 

K  -  PP'I^  ‘■‘P  [■  -  ♦')  ■  'k  PP  ] 

On  using  (14)  it  is  easily  established  fiom  (1),  (13),  (17),  and  (19)  that  asynqitotically  as 

r’/r  oo 


E(p’,z')  ~  ^f(8’) 


(20) 


where  cosO'^z'/r* ,  and  the  "scattering  amplitude"  is 


[exp(-  ikzcosO')  J  Ji(kpsin0’) 


(21) 


In  die  integrand  of  (19)  one  must  evaluate  all  the  derivatives  and  dim  express  p  and  z  in 
terms  of  s. 

3>The  Tapered  Wave^de,  Horn,  and  Scatterer 

The  configurations  with  wl^h  we  shall  be  ccmconed  are  ttqiered  waveguides ,  cavities 
(both  closed  ^  open),  horns,  and  scatterers.Considff  first  a  tapered  perfectly  conducting 
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wav^uide  as  shown  sdianarically  in  fig  1 .  The  metallic  wall  indicated  by  the  heavy  line  is 
assta^  to  have  an  assigned  thickness  and  sh^w.  The  sections  on  die  Car  right  and  left  are 
suf^xised  to  be  long  straight  sections  of  ctmstant  radius  R+ on  the  rig^t  and  ontheleft. 

They  are  joined  by  a  time^  section  die  radius  of  which  is  R  «  R(z).  The  system  is  taken 
to  be  enctosed  mettuHc  walls  winch  are  assumed  to  be  ideally  conducting  This  requires 

diat  the  tangential  cooqxMient  of  die  electrical  fidd  vanish  on  dre  walls,  whrace,  foUosring 

(1),  YaO  on  die  walls. 

For  an  infinite  cylindrical  wavegmde  radius  Rt ,  as  can  be  readily  determined  by 
separatitm  of  variables,  eq.  (15)  has  eigensdutions 


where  n*  1,2,3,- 


A 


2  _ 


JlO'iii)*0 


Suppose  that  ^  is  positive  for  na0,l,2,  —  ,  and  negative  otherwise.  Then 

all  modes  with  n>N±  will  decay  as  one  advances  out  of  the  tqio’ towards  the 
^*j||**^  dying  out  in  a  distance  of  the  order  of  .  Thus  far  to  the  left  of  the 

tapor,  Y  can  be  well  rqiresenied  by  the  finite  series 
N-  N- 

T(p^)  -  y  «;F.(|:]«p(ik;  .y  y  p;F.[£-]«pr-iic;  z)  (26) 

Assume  that  the  source  rtf  the  ratfiation  is  on  the  Idt  This  corre^ionds  to  stipulating  on 
die  segment  F- of  die  boundary  curve  r  diat  is  at  a  point  z- sufficiently  ftff  to  die  left  of  the 
taper  so  diat  the  evanescent  modes  genoated  by  the  taper  have  died  out,  the  coefficients 


7 


of  tberightwaidpropagatiiig  modes.  Moreover  one  mast  require  on  the  segment  Ef  of 

the  boundary  r,  that  is  at  a  fxwt  Zf  sufficiently  far  to  die  ri^t  that  all  the  evanescent 
modes  generated  by  the  taper  have  died  out ,  diat  the  coefficients  ^  of  die  leftward 

pn^iagating  solution  vanish,  corresponding  to  the  absence  of  a  source  on  the  right,  namely 
diat 


N+ 

'l'(p,z)  =  Fn^^jexp[ik„(z-z^)]  (27) 

n=0 

Also  Y  must  vanish  on  the  segment  Fw  corresponding  to  the  perfecdy  conducting  wall  The 
task  is  thus  to  find  die  coefficients  at  and  .  On  the  perfecdy  conducting  wall,  the 

segment  Fw,  one  must  require  that  ^=0.  As  we  shall  see  in  section  4,  this  is  accomplished 
by  decomposing  die  problem  into  a  number  of  simpler  inoblems  such  that  die  answer  is  die 
sum  of  the  solutions  of  these  simpler  cases. 

The  case  of  a  miciowave  hom  is  shown  sctematically  in  Hgure  2.  The  madiematical 

problem  is  the  same  as  that  of  the  tapered  waveguide  except  that  the  curve  F+  has  moved 
off  to  infinity  where  'F  must  vanish ,  which  feature  is  automatic  because  of  the  asynqitotic 

properties  q€  the  Green's  function.  On  the  outside  of  the  hom  the  segment  of  F  there  is 
terrninatBdsufficiendy  far  to  the  left  that  the  diffiacted  field  is  negligible.  As  in  the  case  of 

the  tapered  waveguide  one  has  to  stipulate  the  0^  and  determine  the  .  On  the  perfecdy 

conducting  wall,  the  segment  Fwf  cme  must  require  that  Y=0. 

For  the  case  of  a  scatterer  as  shown  schematically  in  Figure  3,  the  domain  of  interest  is 
that  external  to  the  object  One  must  give  the  incidoit  wave  'Fin,  for  example  a  plane  wave 
e**™-  It  is  then  convenient  to  write  'F=  + 'Fs .  On  the  scatterer,  the  segment  Fw,  where 

die  electric  fidd  and  hence  Y  must  vanish,  die  scattered  wave  'Fj  satisfies 
'5's(P(s).z(s))  =  -  'Fin,  where  s  denotes  arc  length  along  the  curve  F,  the  intersection  of 

the  outside  of  the  scatterer  and  a  plane  of  constant  azimuth  ^.Clearly  one  can  deal  with  the 
case  where  the  scatterer  is  inside  i.  waveguide  or  hom  in  a  parallel  way  by  combining  the 
problems  already  described. 

For  a  closed  cavity  die  domain  is  the  interior  of  die  object,  the  hatched  region  of  Figure 
2.  In  this  case  the  frequency  (d  is  not  determined  by  an  external  agency,  rather  it  is  found  to 
be  an  dgenvalue,  of  which  there  are  an  infinite  denumerable  se^.  It  is  is  readily  seen  that  the 
eigenvalues  oi^  are  all  real  and  positive.  One  multiplies  (5)  by  'F'*  and  integrates  over  the 


volume  of  the  cavity,  whence  on  using  Gauss's  theorem  and  the  requirement  diat  4*  vanish 
on  the  metallic  wall  it  follows,  even  when  e  and  ^  are  functions  of  positicm,  that 


For  the  cavity  the  term  ^s')  in  (15)  vanishes,  the  integral  equation  is  homogeneous,  and 

has  solutions  only  for  characteristic  values  of  oo?. 

For  the  case  of  an  open  cavity,  like  the  opoi  resonator  shown  schematically  in  Figure 
4,  the  associated  inte^  equation  is  homogroeous,  but  die  eigenvalues  are  complex 
bwause  of  radiation  loss.  The  most  interesting  case  of  low  systems  where  foe 
eigenvalues  are  almost  teal,  can  be  dealt  with  very  efficiently  uod  die  Irest  lossy 
eigenvalues  found  efBciendy. 

4.  Reduction  of  Problems  I  and  II  to  Sequences  of  Dirichlet  Problems 
While  foe  {uoblems  presented  in  section  3  are  weU-posed,  foe  analytical  and  numerical 
apparatus  for  dealing  wifo  them  is  not  as  well  devek^ied  as  that  for  Clascal  problems  of 
mathematical  physics,  such  as  Dirichlet  and  Neumann  problems  for  the  Fklnoholtz 
equation^.  Therefore,  we  will  reduce  them  to  finite  sequences  of  Dirichlet  problems  for  (5) 
<Mi  perspicuously  constructed  regions,  and  later  deal  wifo  them  numerically  . 

Consider  the  case  of  foe  hom  depicted  in  Figure  2.  Recall  that  one  stipulates  the  waves 

incident  from  the  left  on  F-  (at  zsz.)  by  giving  foe  coefficients  a~  in  (26),  and  requires 

that  Y  vanish  tm  foe  metallic  wall  corresponding  to  the  curve  Fw  In  ender  to  deal  wifo  this 
situation  and  determine  the  unknown  coefBcients  ^  characterizing  die  reflected  waves  we 

construct  a  finite  sequence  of  auxiliary  functions  4>i  defined  in  foe  dennain  Q  in  p,z  interior 
to  the  curve  FaF-’i-F^  +  the  sector  of  a  ciicfe  the  radius  of  which  tends  to  infinity,  as 
indicated  in  Hgure  2.  The  desired  solution  *P  can  be  exfaessed  as  a  linear  combination  of 

these  and  die  coefiBdents  determined  in  a  convenient  manner.The  <t>j  are  solutions  of  foe 
following  Dirichlet  problem  (where  for  the  hom  {xoblem  we  suppress  any  unnecessary 

superscripts  and  subscripts  -). 

a)  For  all  m  s  i;2,. . .,2N,  tbni  satisfies  Equation  (32)  inside  Q. 

b)  For  all  m  a  1^...^,  tbm  vanishes  on  Fw 

c)  For  all  mal,2,...J^,  onF- 


(29) 
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2N_ 


2N- 


YmAm^  cxp(ik|i)  + 


y.  YmBmji  exp(-ikoz) 


(35) 


ni»l  m=l 

=  On  e*p(iknz>4-  ft,  exp(-iknz) 

The  coefficients  d  expOknz)  and  exp(-ikaz)  in  (35)  must  vanish  separately,  as  can  be  seen 
by  multiplying  (35)  by  exp(±ikaz)  and  int^rating  with  respect  to  z  over  a  distance  2]t/kn- 
Thus  for  n=l,2,3...N_ 


2N- 

m=l 


On 


2N_ 

*  fti 


m«l 


(36) 


Eqs.  (36)  are  a  system  of  2N-  linear  equations  for  die  unknown  N-  coefficients  Ym  and  the 

unknown  N-  coefficients  fta.  With  the  coeffictots  thus  determined  (33)  is  die  desired 
solution  since  it  satisfies  the  differential  equatkm  (6)  and  all  the  boundary  conditions. 

For  the  case  of  the  tapered  waveguide  the  procedure  is  similar.  The  domain  Cl  is  defined  by 
the  z  axis  and  die  curve  rsF-  u  u  (see  Figure  1).  One  defines  auxiliary  functions 
‘I»m  as  follows. 


a)  For  all  m  «  1,2,. .  .,2N-,  ^  satisfies  Equation  (32)  inside  O. 

b)  Forallms  i^„..^N-,  Om  vanishes  on  Fw-t-F.,..  (37) 

c)  For  all  m  a  1,2,...  onF- 


(38) 


d)  For  aU  m  »  N-+1 J4-+2,.  ..,2N..  on  F_ 


exp^-ik^zj 


(39) 


One  defines  auxihary  functions  functions  'Fm  such  tiiat 

^  ^rall  m®  1,2,...^N+,  'Flu  satisfies  Equation  (32)  inside  £2. 

b)  Fot  all  m  =  1^,. .  .^N+,  ^ni  vanishes  on  Tw  +  F— 

c)  ForaU  m=  onr+ 


d)  For  all  m  ®  N++l,N++2,. . .^N+.  on  r+ 


‘I'm=  Fm(|-]  exp^'-ikj^zj 


Then  far  to  the  right  of  the  taper 
N+ 


N+  N+ 


N+ 


N. 


n®l  n=l 


and  far  to  the  left  of  the  taper 


N. 


N- 


n-1  n-1  ^ 


N.  N_ 

'I'mCp.z)  -  ^  ^  m 

n-1  n-1  ^  ^ 

The  coefficients  in  (43)  -  (46)  are  determined  by  evaluating  die  solutions  of  the  integral 

equations  for  <I>ii  and  M  an  appropriate  number  of  suitably  chosen  points. 

Next  one  writes 

2N-  2N+ 

^(P,z)=  Xllin‘*Vn+  (47) 

m=l  m=l 

On  equating  the  asymptotic  repiesMtations  (26)  and  (27)  on  their  respective  domains  of 
validity ,  and  using  the  orthogonality  of  the  Fn  and  integration  over  z  as  was  done  fw  the 
h(xn ,  one  obtains  the  algebraic  system 


2N_ 

2N+ 

(48) 

m=l 

msl 

2N- 

2N+ 

* 

■  X  *  a”(given) 

m,n  n 

(49) 

m=l 

m=l 

2N_ 

2N+ 

YB--+Yn. 

Zmd  ni,n' 

(50) 

msl 

m=l 

2N_ 

2N+ 

X®  Ym  +  B'*”“T|in  =p 

m^'  m,n  n 

m=l  in=l 


(51) 
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Eqs.  (48)  -  ( 51)  are  a  set  of  2(N-+  N+)  linear  equations  for  the  N-  conq)lex  numbers  Ym, 

the  N+  coiiq>lex  nutObers  1110,  the  N-  cotiq)lex  nundrers  pT,  and  the  complex  numbers 

n 

ttg  .  Once  these  have  been  found  the  solution  ^(p^)  is  ccmipletely  known  via  (47),  and 
transmission  and  reflection  are  deteimined  by  tiie  coefficients  0^  and^. 


S'Nystrom  method  for  the  solution  of  second  kind  integral  equations 
(SKIEs) 

We  have  reduced  the  problem  of  determining  die  electric  field  to  consideration  of  (IS), 
a  second  kind  integral  equation  (SKIE).  This  can  be  dealt  with  numerically  by 
discretization  via  an  appiopiiatB  quadrature  focmu^. 

We  will  define  an  n-p^t  quadrature  rule  on  the  interval  [OJ^]  as  a  finite  sequence  of 

pairs  {xi,tt>i),  isl,2,....n,  where  xi€[04<]  fn^all  i  6[l4i].  For  a  function  y : 
[0J^]-^Rl,  we  will  look  upon  the  sum 


n 

n(¥)  =  ]L«>i¥(xi)  (52) 

i=l 


as  an  approximation  to  the  integral 

L 

J  y(x)dx  .  (53) 

For  the  cases  considered,  we  treat  L  is  the  length  of  the  arc  F.  We  will  say  diat  the  family 

of  quadrature  formulae{Tin=  {x^’^n  )  ’  *  *  lA-.^i} ,  n  =  1,2,...,  has  a  convergence 

rate  m  (m^l)  for  the  function  y :  [0J.r]--»Rl>  if  there  exist  numbers  A>0  and  integers  N>0 
such  that 


L 

I  VV)-Jv(x)dxl<^ 


(54) 


for  all  n  >  N. 

In  order  to  solve  the  integral  equation 

L 

o(x)  +  j  F(t,x)  a(t)  dt  =  <|>(x)  (55) 


the  Nystrom  algorithm  tqilaces  (32)  with  a  system  of  linear  algebraic  equations 


(56) 


n 

y  F(xj^  =  ♦(xi), 

with  i  =  1,2,... 41,  and  10^,  to^,  ...,to°tbeooe£BdeDtsofanai]pioptiatd[y  chosen 

quadrature  fonnulae.  We  win  denote  the  matrix  of  the  system  (27)  by  An,  and  view  die 
solution  Oi,(^,...,an  of  (27)  as  an  approximation  to  die  solution  ^9  of  (26)  at  die  nodes 

12  n 

Note  that  in  general,  the  ooeffidentstOj,  t0j,...,T0°  of  the  quadrature  jfbrmula  depend 

on  the  pdnt  Xi  at  which  the  integral  (32)  is  bemg  qipioximated.  When  the  kernel  F  is 
smooth,  this  can  be  avoided  by  choosing  a  sing^  quadrature  formula,  such  as  the  end-point 
corrected  trapezoidal  rule,  that  win  perform  u^iwaU  is  1,2,.. .40.  However,  in  many 
applicadons  the  kernel  is  singular,  and  die  distribution  of  its  singularides  usuaUy  fmres  a 
qiecial  choice  of  aquadrature  formula  for  each  node. 

The  fdlowing  theorem^  is  the  princ^  justification  for  the  use  of  the  Nystrom 
algorithm  for  die  solution  of  SKIEs .  Siqipose  diat  Equation  (32)  has  a  unique  solutirm, 

and  that  its  kernel  F  and  right-hand  side  ^  are  continuous.  Suppose  futther  that  a  family 
(iln)  of  quadrature  formulae  is  such  diat  for  some  B  >  0 


n 


I 

isl 


(57) 


forallnsi^,,,. .  Then  diere  exists  a  numb^OO  such  that  for  any  n^Zand 

i  *  ltZy.»»tnt 


(58) 


where  k  is  the  convergence  rate  of  the  fonnulae  {tIq}  . 


5.  Dcscriptioii  of  the  Algorithms 

Now,  armed  with  die  qiparatus  developed  in  the  preceding  sections,  we  are  prepared  to 
construct  algorithms  for  die  solution  of  the  original  proUems  in  Section  3.  Fir^we 
describe  an  algorithm  for  the  solution  of  the  htm  problem. 
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Algorithm  I 

Step  1.  Construct  the  curve  F  by  appending  the  segment  F-  to  the  user-specified 
waveguide  boundary  (see  Figure  5).  Use  spline  package  FTTPACK^  to  construct  an  equi- 
spaced  discretization  of  the  resulting  curve.  Select  the  N  nodes  into  which  the  boundary  F 
is  to  be  discretized. 

Step  2.  Convert  the  problem  into  a  sequence  of  exterior  Dirichlet  problems  in  the 
region  n  as  described  in  section  4.  F<x  each  i  »  1^,...,2N,  convert  the  corresponding 

Dirichlet  problem  into  a  second  kind  integral  equation  (SKIE)  on  F. 

Step  4.  For  each  i  «  1,2,.. .,2N.  discretize  the  obtained  SKIE  via  the  Nystrom 
method  based  on  the  quadrature  formulae  of  Appendix  C,  obtaining  a  system  of  linear 
algebraic  equations  of  (umension  2N  . 

Step  5.  For  each  i «  1,2,.  ..,2N,  solve  the  linear  system  obtained  on  the  preceding 
step  by  means  of  a  standard  Gaussian  eliminadoi  subroutine  firom  UNPACK  ^ .  View  the 
solutitm  as  an  apixoximation  to  o  on  F  solving  the  underlying  Dirichlet  problem. 

Step  6.  Find  the  coefficients  Yi.Y2.'">Y2N  in  equation  (33)  by  solving  (36) 

Step  7.  Linearly  combine  the  soluti<»s  of  the  2N  linear  systems  obtained  on  the 
preceding  step  according  to  the  fcvmula  (33).  The  result  is  a  distribution  o  on  Y2  whose 
field  is  a  solution  of  the  problem . 

Step  8.  For  each  point  x€R^  where  the  solution  of  the  horn  problem  is  desired, 
evaluate  it  by  approximating  the  integral 

J  o(t)  G(x,t)  dt  (59) 

F 

by  the  trapezoidal  sum 
N 

h  I  o(ti)  G(ti,x)  (60) 

i=l 

with  {ti}4  =  1,2,  —  N,  the  nodes  in  the  discretization  of  F,  and  h  the  sampling  distance 
between  the  adjacent  no^  q,q.i. 

Next  we  dncribe  an  algorithm  for  the  solution  of  the  tapered  waveguide  problem . 
Algorithm  2 

Step  1.  Construct  the  curve  F  by  tqjpending  die  segments  F+,  and  F_  to  the  desired  part 
of  the  boundary  of  the  usn-specified  waveguide  F  (sec  Figure  1).  Use  spline  package 
FTTPACK  ^  to  construct  an  equispaced  discretization  of  the  resulting  curve.  Select  the 
number  N  of  nodes  into  which  the  boundary  F  is  to  be  discretized. 

Step  2.  Convert  the  problem  into  a  pair  of  sequences  of  interior  Dirichlet  problems  in 

the  region  D. 


Step  3.  Foreach  1^,...^-,  convert  die  corresponding  Dtrichlet  Problem  intoa 
second  kind  integral  equation  on  F. 

Step  4.  Foreach  is  1.2....^^  convert  die  corre^)onding  Diiichlet  Problem  intoa 
second  land  mi^nl  equation  on  F. 

Steps.  For  eadi  {•  discretize  die  obtained  SKIE  via  die  Nystiom 

method  based  on  die  quadra^  fonnnlaeofAppendixC,  obtaining  a  systnn  of  linear 
algehndcequatioos  of  dimension  N 

Stqpd.  Foreach  i*  1,2,....2N-),  discretize  die  obtained  SKIE  via  the  Nystrom 
method  bued  on  the  miadratnie  fixmt^  of  Appendix  C,  obtaining  a  system  (rf  linear 
algdxaic  equadons  of  dimension  N. 

Step?.  Foreach  i*  l,2,....2N-,s(dve  the  linear  system  obtained  on  the  preceding 
step  by  means  of  a  standard  Gaussian  elimination  subroutine  from  UNPACK  View  the 

solution  as  an  tq>proKimatk»  to  the  density  a  on  F  solving  die  underlying  Diiichlet  PioUan 

Steps.  Foreach  is  1^...,2N1k  solve  the  linear  system  obtained  cm  the  preceding 
step  by  means  (tf  a  standard  Gaussian  elimination  subroutine  from  UNPACK.  View  die 
solution  as  an  approximation  to  die  density  of  dipole  disttibudon  on  F  stdvingthe 
underlying  Dirkhlet  proUem. 

Stq)9.  HndtheooefficientsYi,1i;2,”‘*‘f2N-*’*^^l.%‘'*»^2i^- 

SiqilO.  Unearly  combine  the  sedations  of  the  2(N--fN4.)  linear  system  obtained  on 
the  preceding  stq;>  according  to  the  ftxmula  (47).  The  result  is  a  distribution  o  on  F  whose 
field  is  a  solution  of  the  tapmd  waveguide  problem. 

Step  1 1.  For  each  point  x€Q  where  die  solution  of  problem  is  desired,  evaluate  it  by 
approxhnating  die  inte^ 

/  o(t)  F(x.t)  dt  (61) 

F 

by  the  trapezoidal  sum 

N 

h  o(ti)  F(ti.x)  (62) 

with  (q),  i  s  1,2,.. .,N  the  nodes  in  the  discretization  cf  G3,  and  h  the  san^ling  distance 
between  the  adjacent  nodes  q,q.i. 

ti’Results 

A  fortran  program  has  beat  written  inqdementing  algorithms  I  and  n  of  this  paper.  The 
program  has  bem  teased  on  a  variety  of  protdems.  We  present  a  detailed  descr^on  and  die 
remits  of  four  such  numerical  experiments. 


Exainirte  1  -  Tqiered  Waveguide 

The  waveguide  defined  in  fig.  5  has  a  left  hand  asynqjtodc  radius  of  0J24  cm.  and  a 
ri^t  hand  asyn^piotic  radius  of  0.743  cm .  Ihe  slu^  was  assigned  stipulating  a  dense 

set  of  points,  wmdi  wete  dim  fitted  to  a  ^line,  a  convenient  dragn  procedure.  The 
transinon  distance  is  roug^y  US  cm.  The  operating  frequency  chmen  is  3.2Xl0il  sec*^ 
This  leads  to  one  propanting  mode  on  die  m  left  ra  the  tape^  section ,  and  two 
propagating  modes  on  me  far  ri^L  The  boundary  enqih^^  for  the  calculation  was  first 

approximated  as  shown  in  fig.  6,  adding  the  pecewise  linear  curves  and  F- to  provide  a 

boundi^'inoo^  die  z  axis.  This  shape  was  dien  resampled  to  pcov^  the  axioodi  curve 
fig.7,  which  was  then  diaaetized  into  90  rouf^y  equally  spaoM  nodes.  Also  indicated 
are  tte  test  points  enqilqyed  in  confuting  the  modal  eiqiansioo  of  die  solution.  The  results 

fra*  the  lines  of  constant  real  part  of  4*  are  di^layed  in  fig.  8.  The  calculation  todc  six 
minutes  of  CPU  time  on  a  VAX  780. 

Exanqile  2  -  0^  Wav^uide 

T^  calculation  usesasemi-infinite  waveguide  of  2  cm.  inner  radius  and  0.25  cm. 

duckness.  The  operating  frequency  is  tixlO^^  sec~^  resulting  in  one  propagating  mode  on 
the  fin  left.  The  lesaaqMedboundi^  is  shown  in  fig.  9.  This  was  discrttized  into  180 
nodm,  with  the  points  most  dense  in  the  neighborii^  of  the  return  in  the  boundary  curve 
at  the  end  of  the  waveguide .  The  calculation  took  17  minutes  of  CPU  time  on  a  Vax  8600. 

The  resulting  Hnes  of  constant  real  part  of  Y  are  shown  in  fig.  10.  The  resulting  radiation 

pattern  is  shown  in  fig.  (1 1)  where  If(6)i2  a  IE|2fr2  is  plotted  vs.  the  pcdar  angle  6  diat  die 
direction  of  observation  makes  with  die  axis  if  the  wav^uide. 

Example  3  -  Microwave  Horn 

Tte  internal  fiur  left  radius  of  the  Ikho  is  2  cm.,  and  its  waU  thickness  0.25  cm.  The 

frequency  enqiloyed  is  6x10^9  sec~i,  yielding  one  propagating  mode  on  die  fin  left.  The 
resa^M  boundary  is  shown  in  fig.  12.  It  was  discretu^  into  180  nodes,  with  the  nodal 
doisity  greatest  near  the  return  in  the  curve.  The  calculation  took  17  minutes  of  CPU  time 

(Ml  a  Vax  8600.  The  lines  of  constant  imaginaiy  part  of  Y  are  shown  in  fig.  13.  The 
radiation  pattern  is  displayed  in  fig.  14,  where  lf(6)P  ^  IEP/r2  is  plotted  vs.  the  polar  angle 
6  that  the  direction  of  obervation  makes  with  the  axis  if  the  waveguide. 

Examine  4  -  hijection  TYough  a  Ifole  in  a  Metallic  Wall 

Tm  fin  left  radius  of  the  waveguide  is  3  cm.  and  its  wall  thickness  is  2  cm.  The 

frequency  employ  is  6x1019  sec~i,  yielding  one  proi>agating  mode  (mi  the  far  left.  The 
resampM  bouiKuny  is  shown  in  fig.  17.  It  was  di«netized  into  180  nodes.  The  calculation 

took  13  minotes  on  a  Vax  8600.  The  lines  of  constant  real  part  of  Y  are  shown  in  fig.  17 

and  die  lines  of  constant  imaginary  pan  (rf 'Fare  (delayed  in  fig.  17.  The  radiation  pattern 

isdi^li^edinfig.  18,  where  lf(6)PslEP/r2  is  {dotted  vs.  the  po*ar  angle  6  that  the 
direction  of  observation  makes  widi  the  axis  if  the  wavi^uide. 

Appendix  A  •  The  Green's  Function 


Recall  that  R  *  Ir  -r1  *  [p2+p'^2pp'cos(^-4')+(*~*')^]^  where  p,^,z  are 
cylindrical  coofdinates.  Then  the  three  dimension  Green's  function 


satisfies 


n 


V2G+k2G  =  -  8(r  -  O  -  JKp.-p:)g(»-  .z’)8(»-»’) 

P’ 


(Al) 


(A2) 


Write 


G(p.p'.<|H»’^-x')  =  F(p.p'.^'.z-a')  €»(♦-♦•) 


Then  it  is  readily  calculated  on  uring  the  propenies  of  die  I^rac  delta  function  that 


(A3) 


S(p  -  p')8(z  -  z')8(»  -  »■) 


Let 


2k 

F(p,p’,z-z’)  =  Jd^  F(p,p’,^^',z-z’) 


Then  it  follows  from  (A4)  on  integrating  over  one  period  in  ^  that 

r-?y 


+  fk2  _  1.>|F  - 

p’ 


Define 


2x 


(A5) 


(A6) 


K(p,p',z-z')  *  pp’  F(p,p',z-z')  *  pp’Jd^  e->(^-^')G(p,p',^',z-z')  (A7) 

Note  dut  K  is  invariant  under  interchange  of  primed  and  unprimed  arguments.  Then  (AS) 
implies 


(A8) 
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But 


^  V^K  2  aK 

(p2  J"  p2  p3ap 

while 

V2r^v— -~+- 

\pj  p  p2ap  p2 


whence  (A8)  in^lies 


pV-  ~  +k2i^=  -8(p-p’)8(z-z*) 


and  on  using  the  properties  (tf  die  Dirac  delta  fonctioa 


L{K)  - 


p2  I  p2 


S(p-p-)6(z-z') 

(pp')»/2 


(A9) 


(AlO) 


(All) 


(A12) 


We  wish  to  derive  an  integral  rcinesenQUion  for  Y  the  fonn 

«P' .  >F(p',rO  =  Jds  CW  [n.VK(p,p'.r-z')]p_p^^^  (AI3) 

r 

where  F  is  that  curve  which  when  rotated  around  the  z  axis  generates  the  surface  S 
enclodng  die  volume  V  containing  the  point  f.  The  element  of  arc  length  along  C  is  ds. 
Note  that  in  virtue  of  its  invariance  under  interchange  of  primed  and  unprimed  arguments  it 

fdlows  from  (All)  that  K  satisfies  L'(K}sO,  wtence  it  is  evident  that  'V  as  defined  by 
(A16)  satisfies  L'{H''}  at  all  points  p',z'  interior  to  the  domain  enclosed  by  F. 

In  (xder  to  derive  an  integral  equation  for  the  function  o,  assuming  diat  Y  is  given  on 
F,  we  consider  a  point  on  F  corresponding  to  s',  and  a  point  p',z'  inside  die  area  in 
question  a  distance  e  from  F  as  measured  akmg  die  unit  outward  ptmting  normal  n(s')  as 
shown  in  Figure  19 .  Hck  e  «  S  «  die  radius  of  curvature  at  s',  and  break  the  domain  of 
integration  up  into  a  segment  of  width  28  straddling  die  pdnt  s',  which  segment  may  be 
cmiadeted  locally  flat,  and  the  rest  of  F. 
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As  p',z'  approaches  F  one  sees  ftxxn  (A7)  tfiat  the  integrand  divei^ges.  The  major 
contributions  to  die  integral  come  fiom  the  deotxninatar  in  the  neighborhood  of  where 
R2»(p-p')2  On  keeping  only  dominant  teims  it  follows  diat 


Thus  the  contribution  to  (A15)  of  the  segment  of  width  26  can  be  written  to  good 
i^ipioximation  as 


s'+8 

s' 

g 

=  2o(s’')p'aictan- 
e 

Qeariy  in  the  limit  e-K)  one  has 


I»  2)kj(s'’)p* 


(A16) 


Thus  in  the  limit  that  the  point  p',z'  lies  on  C  eq.(lS)  yields 


2y(p(s’) 


=  o(s’)  +  P  Jds  <j(s)  |n.V  I 


p*p(s)  z=z(s) 


(A17) 


where  P  denotes  the  Cauchy  principle  value  of  the  integral.  On  introducing  the  definitions 
of  (21)  and  (23)  this  is  die  desired  integral  equation. 
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Fitpireraprions 


Fig.  1  Sdiemadc  diagram  of  a  tiq)eml  waveguide. 

Bg.  2  SdiematicdiagnmofalionL 

Fig.  3  Schematic  diapam  of  a  cavi^  or  scatterer.  For  the  scattercr  the  domain  is 

the  exterior  of  the  boundary  curve,  for  the  cavity  die  domain  is  the  cross 
hatdied  interior. 

Bg.  4  Schematic  diagram  of  an  r^ien  resonator. 

Bg.  S  The  waveguide  example  1  as  ^ledfled  by  the  user. 

Bg.  6  The  wav^uide  of  example  1  after  ixeHminarygeQCDMricprocesring. 

Bg.  7  Resanqiled  waveguide  with  test  points  for  exfuqik  1. 

Bg.  8  EquqioMtials  Of  die  teal  part  of  foe  solution  of  example  1. 

Bg.  9  Resainided  waveguide  with  test  points  for  exan^le  2. 

Bg.  10  Equmotentials  of  the  imaginary  part  <rf  the  solution  to  example  2. 

Fig.  11  Far  field  anplitude  as  a  function  (Wangle  for  example  2. 

Bg.  12  Resampled  waveguute  and  test  points  ftv  example  3. 

Fig.  13  Equmotendals  for  the  imaginary  part  of  the  solution  for  exangle  3. 

Fig.  14  Far  field  anqriitude  as  a  filiation  of  ange  for  example  3. 

Bg.  15  Resaoqiled  waveguide  and  test  pmnts  for  example  4. 

Bg.  16  Equipotentials  for  foe  real  part  of  the  solution  for  example  4. 

Fig.  17  Eqogotentials  for  the  imaginary  part  of  the  solution  for  example  4. 

Bg.  18  Fi^  anglitude  as  a  function  of  angte  for  exangle  4. 
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Figure  6 

The  Waveguide  in  Example  1  After  Preliminary  Geometrical  Processln 


Figure  7 

Resampled  Waveguide  with  Test-  Polnis  In  Example  1 


Figure  8 

Equlpotenilal  Lines  of  the  Reel  Port  of  the  Solution 

in  Exofflple  1 


Figure  11 

Far-Fleld  Pmplliude  os  a  Function  of  Angle 
in  Example  2 
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Figure  12 

Resampled  Waveguide  with  Test  Points  in  Example  3 


Figure  13 

Equlpotent lal  Lines  of  the  Imaginary  Part  of  the  Solution 

In  Example  3 
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Figure  14 

flmplliude  as  a  Function  of  Angle 
in  Example  3 


Figure  15 

Resampled  Waveguide  with  Test  Points  In  Example  4 
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Figure  16 

Equlpotent lal  Lines  of  the  Real  Part  of  the  Solution 

In  Example  4 
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Figure  17 

Equlpotent lal  Lines  of  the  Imaginary  Part  of  the  Solution 

In  Example  4 
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Figure  18 

Far-Fleld  Amplitude  as  a  Function  of  Angle 
In  Example  4 


